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Abstract 

We propose a new method for an analytical, non-perturbative computation of ef- 
fective quark interactions from QCD. It is based on an exact flow equation which 
describes the scale dependence of the effective average action for quarks in presence 
of gluons. 



1 Introduction 



Quantum chromodynamics as the theory of strong interactions shows very different 
facets at short and long distances. The high momentum behaviour is governed by 
asymptotic freedom |I| and the relevant degrees of freedom are quarks and gluons. 
In contrast, the particles which are observed at large length scales (larger than 
1 fm) are mesons and hadrons. The interactions of the pseudoscalar mesons are 
modeled by chiral perturbation theory 0. The corresponding nonlinear a-model 
shares the flavour symmetries of perturbative QCD. It is believed that the free 
phenomenological parameters of this model can ultimately be computed from the 
action of QCD, but this is not a simple task. 

Recently, the transition from quark degrees of freedom to meson degrees of free- 
dom has been described by a nonperturbative flow equation ||^. It is based on the 
concept of the effective average action Tk [0 for which only quantum fluctuations 
with (covariant) momenta > k"^ are integrated out. The average action is the ef- 
fective action for averages of fields. It acts like a microscope by which we can look at 
the theory at different length scales, with a "resolution" given by the scale k~^. For 
k = one recovers the usual effective action, i.e. the generating functional of the IPI 
Green functions, whereas for k —>■ (yo equals the classical action. The dependence 
of Ffc on the scale k is described by an exact nonperturbative flow equation |p. By 
the introduction of composite fields for the mesons the original equation formulated 
for quarks can be transmuted into an equivalent exact flow equation involving also 
the mesons 0. 

In a first attempt the chiral condensate < ipip > and the pion decay constant 
f-jr were computed along these lines 0. In this approach the gluons have not been 
considered explicitly. Their effect was encoded in a phenomenologically motivated 
four-quark interaction, which may be thought as the result of integrating out the 
gluons in the defining functional integral. If one aims at a computation of the effec- 
tive parameters of chiral perturbation theory from the QCD action, this shortcut has 
to be removed. A full computation should start at short distances with the quark- 
gluon description of QCD and systematically account for the quantum fluctuations 
of the gluon field. 

Conceptually, one may integrate out the gluons at once and end with an effective 
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quark theory. In practice, this seems almost impossible: A determination of the 
resulting complicated nonlocal quark interactions requires more or less a complete 
solution of QCD. On the other hand, the form of the flow equation at a given scale 
k is not sensitive to details of physics at momentum scales much below k. Only 
the effective propagators and vertices for modes with momenta ~ k'^ play a 
role for this equation. A typical scale where the light mesons form is aroimd 700 
MeV and therefore substantially higher than the confinement scale. This situation 
further improves for problems involving heavier quarks as for example the charmed 
or beauty mesons. One may hope that very detailed features of confinement are not 
needed for an understanding of the mesons. An appropriate tool would therefore be 
a method which only integrates out the gluon fluctuations with momenta > k^ 
instead of addressing the much more complicated problem of integrating out all gluon 
fluctuations at once. Formally, it is easy to do this at a given scale ki. One needs to 
compute the effective action for quarks and gluons FfeJ'^, A] at this scale. Solving 
the classical fleld equation for the gluon fleld A in dependence on ip and inserting 
this classical solution into Ffeif'^, A] yields exactly an effective average action Fj^J'^] 
involving only the quark fields. The quark-effective action Fj^^ [■0] may then be used 
as an initial value for solving a pure fermionic fiow equation for k < ki. Obviously, 
the shortcoming of such an approach ist the complete omission of the effects of 
gluon fluctuations with < kj. Choosing a different scale ^2 for eliminating the 
gluons will lead to a different result. Such a sharp transition between the quark- 
gluon system and a description involving only quarks necessarily introduces a certain 
degree of arbritariness. 

We propose here a more refined method which changes the classical solution for 
the gluon field in the course of the evolution towards lower k, thus refiecting the 
change in the form of F^ [■(/', A]. The result is a smooth procedure for integrating out 
the gluons, where at every scale k all gluon fiuctuations with q"^ > k'^ are included. 
Nevertheless, we obtain a fiow equation for the quark effective average action FfefV'], 
where gluon fields do not appear explicitly. As a consequence of the inclusion of 
contributions from additional gluon fiuctuations as the scale is lowered, correction 
terms appear in the flow equation for the quark interactions. In particular, we choose 
here a formulation where the /c-dependent classical solution for A as a functional of 
■0 includes an effective infrared cutoff ~ /c. As a consequence, the only nonlocalities 
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in Tk[ip] concern length scales shorter than k^^. For the fermionic low momentum 
modes [■?/;] is an effectively local action. For example, a derivative expansion 
is meaningful for <^ k'^. The expected nonlocalities arising from the complete 
elimination of gluon fields (for example a four-quark interaction ~ ^ appearing 
already in the Born approximation) build up only step by step as k is lowered to 
zero. As a result of this method we will end with an exact nonperturbative flow 
equation for the scale dependence of F^ [■?/']. Approximations will be needed to solve 
this equation but they are not limited to perturbative concepts. The correction 
terms reflecting the gluon fluctuations require limited knowledge about the effective 
gluon propagator and vertices. It is hoped that rather crude approximations for the 
gluonic vertices can already lead to satisfactory results. 

In sect. 2 we first demonstrate our formalism for a simple model of two scalar 
fields. The flow of the effective average action for one of the scalar fields obtains by 
integrating out the other scalar field at any scale k. Subsequently this is general- 
ized to quarks and gluons. We also give a first demonstration how this formalism 
describes the flow of the two- and four-point function in the effective quark theory. 
The special case of heavy quarks is addressed in sect. 3. Here we argue that the 
evolution of the gauge field propagator is needed in this limit. In sect. 4 we com- 
pute the corresponding flow equation and discuss the scale dependence of the gluon 
propagator. In sect. 5 we collect all the ingredients needed for the flow equation in 
case of light quarks. Finally, our conclusions are contained in sect. 6. 

2 Reduction of degrees of freedom 

First we consider for simplicity two types of scalar fields, (/? and ip. We want to 
develop a formalism how to translate evolution equations for the effective average 
action for (p and into corresponding equations involving only ip. The reader may 
associate </? with the gluon fields and ip with the quark fields. Our aim is then the 
construction of the effective average action for quarks out of the coupled quark-gluon 
system. This amounts to integrating out the gluonic degrees of freedom represented 
in the simplified model by (f. We start with the scale-dependent generating func- 
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tional for the connected Green functions 

Wu[J, K] = \nJ D^'Dij'exp -{S[^', + aI^^ S[^'] + A^ S[^P'] - JV " ^V'} 

(2.1) 

Here we denote the degrees of freedom contained in (f' (for example the Fourier 
modes) by v?'" and similar for ip', J and K, withQ JV = J^V^'", K^ip' = K^ip'^ . 
We have introduced an infrared cutoff quadratic in the fields 

^t^SW] = (2.2) 

and similar for ip. This suppresses the contribution of fluctuations with small mo- 



menta q < k to the functional integral ( |2.1| ). Typically R)^ , RjJ^ are functions of 
as, for example, 



l-exp(-£) 

which acts like a mass term rI^^ ~ Z^/t^ for q"^ <^ ^ . The effective average action 
Ffcfv?, ^Z'] is related to the Legendre transform of Wk[J, K] 

ffc [ip,ij] = -Wk [J, K] + Jt(^ + i^t^ (2.4) 
by subtracting the infrared cutoff term 

r,[y,, V'] = f,[y,,^] - A^;f^s[ip] - (2.5) 

For — > the infrared cutoff A^S = A^^^ S + A^jf^ S vanishes and Fq is the usual 
generating function for the IPI Green functions. Using the quadratic form of AkS 
it is straightforward to derive an exact non-perturbative evolution equation for the 
dependence of the effective average action on the scale k (t = \nk) p 



dTi, 1„ I .-(o). 1 dRh 



:2.6) 



Here f^^^ = F^^^ + Rk and the inverse propagator F^.^"* is the second functional 
derivative of F^ with respect to the fields. The matrix Rj. = R^^^ + R^^^ is block 
diagonal in (p and tjj spaces. The presence of the infrared cutoff Rk in F^ guarantees 
infrared finiteness for the momentum integral implied by the trace even in case of 



"'^We use indices a, a' etc. for ip and /3, (3' etc. for 
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massless modes. Ultraviolet finiteness is guaranteed by the exponential decay of 
dRk/dt (|2.3| ). A solution of the flow equation (pl6|) interpolates between the classical 
action for A; oo (or k equal to some ultraviolet cutoff A) and the effective action 
for — i> 0. 

The generating functional for the connected Green functions for ip obtains /^from 
(O for J = 

Wk[K]=Wk[J = 0,K] (2.7) 
Correspondingly, we may introduce an effective action expressed only in terms of ip 

Tk[ij] = Tk[^,[ij],i;] (2.8) 



r4^] = hm-Afs[4^] 

= T,Mi;],i;] + A^^^S[Mi^]] (2.9) 
by inserting the fc-dependent solution of the field equation 

= (2.10) 

This defines ipk as a fc-dependent functional of tp. It is easy to verify that Tk[ip] is 
the Legendre transform of (|2.7| ). One concludes for /c — that Fq ['?/'] is the 

generating functional for the IPI Green functions for tp. 

We want to employ the flow equation ( p.6|) for finding the fc-dependence of rfc[?/^]. 
In addition to the corresponding equation for only one type of fields we have here 
additional contributions from the fc-dependence of A^^^ S in ( |2.1| ). The evolution 
equation for Tk[ip] can now be obtained by performing in eq. (2.9) a variable trans- 
formation which amounts to a shift of <f around (pk[4']y^'^ = ^p'^ — '^'^\ip\. One 
obtains 




It is easy to verify that this equation reduces in the limit I^^^ = to the equivalent 
of eq. (p.6|) for fields ip only. The corrections in the first two terms involve the 
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explicit form of the "classical solution" v^fei"^]- If one is interested in IPI Green 
functions for ip with a given number of external legs one only needs a polynomial 
expansion of ^k\4'\ to a given order. For example, the evolution of the term ~ if)'^ 
in Tk[Tp] needs the classical solution up to the order if the series v^a:[V'] starts 
with a term quadratic in if). Additional knowledge of the form of rfc[(^, beyond 
its value for = is only needed for the last correction term in the form of 

Only the ^-dependence of the effective (^-propagator plays a role for the study of 
IPI functions for ip. 

Let us next apply the general flow equation explicitly to the quark-gluon 



a- F)rn*rlrn^ \^=0 + [Rk ) „ (2.12) 



system. If is a Grassmann variable as appropriate for fermions the matrix Rk in 
( p.6| ) becomes 

Ru = Rt^-Rf^ (2.13) 

Also ip* should be replaced by ip and the index summation over j3 should involve 
both ip and ip separately. For the gauge fields we will choose here a formulation with 
explicit ghost variables in close analogy, but slightly different from the formulation in 
ref. 0. This makes our formulation as close as possible to the language of standard 
perturbation theory. Details can be found in the appendix. 

We start with the action including a gauge-fixing term in the background gauge 
and a corresponding action for the anticommuting ghost fields ^, ^ 

S[^\i\a-A] = S[ij\A'] + S,f[a-A] + Sgu[i\a-A] (2.14) 

Here 5* is a gauge invariant functional of the fermion fields ip, ip and the gauge field 

^^ = A^ + a^. (2.15) 

The background gauge field A^ appears in the gauge fixing and ghost terms 

S,f = ^Jd'xGlG^ (2.16) 
G' = mA]ryal (2.17) 
S,n = Jd'xey{-D>^[A]D,[A + a])y£^ (2.18) 
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Here -D^[^] is the covariant derivative in the adjoint representation in presence of 
the background gauge field A. The generating functional for the connected Green 
functions is defined as usual 

W[7], C, K;A]= J V^'V^'V^'V^'Vaexp -{S - J d'^x[f]%l)' + r?^' + CC' + Cl' + Ka]} 

(2.19) 

We have introduced here also sources C for the ghost fields and note that the source 
couples to the gauge field fiuctuation and therefore transforms homogeneously 
under gauge transformations as an adjoint tensor. The fc-dependent version 
obtains from W by adding to S the infrared cutoff piece 

Afc^ = A?)^ + A^^^S + aI^'^)^ (2.20) 

Here the fermionic cutoff reads 

A?)^ = 4(/2f )V 



d'xi;'Z^^SYD^[A]) r'^\-D'[A]lk'W (2.21) 

with the covariant derivative in the appropriate representation [D"^ = D^D^) 
and r^^^ a dimensionless function. For the gauge field cutoff we choose 



V[A] r 



{A) 
k 



'z~aI,v[A\ 



< (2.22) 



with T>\A\ an appropriate operator generalizing a covariant Laplacian in the adjoint 
representation which will be explained below. The matrix Z^^k accounts for an 
appropriate wave function renormalization. Finally, we take for the ghosts 

Af)5 = e;,(i?f^)^;e'^ 

= / d'^QZ,K^'^s\A\T'i^\VAm^)\\e (2.23) 

with Ps[A] = — D^[y4] in the adjoint representation. A good choice for the dimen- 
sionless function is0 

rk{y) = (2.24) 



^The function r^jf^ may be chosen differently from (2.24) in order to avoid that Rk diverges for 
vanishing covariant momenta. 
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such that 

hm Rk = Zkk"^ (2.25) 

The fc-dependent functions Z^^k,2A,k and Zgh^k will be adapted to corresponding 
wave function renormalization constants in the kinetic terms for the fermions, gauge 
fields and ghosts. In principle, they can depend on the background field A. The 
infrared cutoff piece AkS cuts off all quantum fluctuations with covariant momenta 
smaller than k in the functional integral defining Wk- For covariant momenta larger 
than k the infrared cutoff is ineffective and its contribution to the propagator is 
exponentially suppressed. 

Performing a Legendre transform and subtracting the IR cutoff piece again (c.f. 
(p.4[), (|2.5|)) we arrive at the effective average action T^lip, ^, A, A], where A = A + a 
and a is conjugate to K. The dependence of on the scale k is described by 



an exact evolution equation analogous to eq. ( p.6|) , with a negative sign for the 
contributions ~ R^^^ and -R^^''^ It is derived in the appendix (A. 12). We note that 
Ffc only involves terms with an even number of ghost fields due to the symmetry 
— — > of the Sgh and A^^^^S. In consequence, the ghost field equations 

^ = 0, ^ = (2.26) 

have always the solution ^ = ^ = 0. We therefore can extract the propagators and 
vertices for the physical particles from the effective action for ^ = ^ = 0: 

Tk[^p,AA] = Tk[^p,0,A,A] (2.27) 

Nevertheless, the evolution equation for Tk[ip, A, A] obtains a contribution from the 
variation of the infrared cutoff of the ghost fields as given by 



-Tr jl^^i^r J (rr^ + i?. j ^-Sk (2.28) 
Sk = Tr{g4^^))(rr(n/?fy^} (2.29) 

Here T^^^ +Rk in (|2.28|) is the matrix of second functional derivatives of FfcH- A^f ■'S' + 
A^f^ S with respect to ip and A at fixed A. As compared to the more symmetric form 
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of the flow equation (A. 12) we have combined here similar pieces in the quark and 
ghost sector. One should remember, however, that the fermionic part (r^.^"* + Rk)~^ 
is a submatrix of a larger matrix containing also ipip entries. For the derivation 



of eq. ( |2.28| ) we have exploited that the matrix of second functional derivatives of 
(Ffc + AkS)[ilj,^, A, A] is block diagonal in the (ipyA) and ^ components for = 
0. The ghost dependence of Tk[ip,^,A,A] appears in the evolution equation for 
Tf:[ip, A, A] only through the term which involves the second functional derivative 
with respect to the ghost fields r^f''*-**-^-', which is evaluated at ^ = ,^ = and may 
depend on ip, A, A. As a consequence of local gauge invariance the average action 
must obey anomalous Slavnov- Taylor identities which are displayed in the appendix. 
They constrain, in particular, the ghost dependence of F^,. We will not pay much 
attention to the detailed form of r^^'^'"^^^ in the present paper and approximate it by 
its "classical" value (cf. ( prTSD ) 

^igh){2) ^ (2.30) 

or a slight generalization thereof (cf. eq. (4.8)). In order to complete the formal 
setup of our investigation we need to specify the operator V in eq. (|2.22|) . A good 
choice is 

V[A] = Ti'^^^^\A] (2.31) 

where F^'^'**'^'' is the second functional derivative of Tk[ip, A, A] with respect to A for 
fixed A and = 0, evaluated at the point A = A. As in previous formulations 0, the 
effective average action F^ [■?/', A, A] is gauge invariant with respect to simultaneous 
gauge transformations of ip,A and A. 

We can now apply the formalism of the last section in order to "integrate out" 
the gluon fields A. The classical field equation, whose solution is Ak, reads 



where the derivative should be taken at fixed A. At this point Af^ becomes a func- 
tional of ip and A. For the purpose of the present paper we only consider the special 



choice A = and omit the argument A in the following. In this version rI^^ be- 



comes a simple function of momenta. Summarizing our adaptation of eq. (|2.11|) for 
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quarks and gluons one obtains 






this equation we need Afc[V^], [t/^, A^] and €k[ip, A^]. 

Our aim is the solution of the flow equation for k —>■ 0, starting at some high 
scale ko where T^g [ip] can be reliably computed by solving the field equation for 
in a perturbative context. As k decreases, we gradually explore the non-perturbative 
regime and the full quantum-effective action obtains for A; = 0. Obviously, such a 
program is only feasible with approximations that truncate the most general form 
of Ffc. The truncation used in the next two sections concerns the three-gluon and 
four-gluon vertices which are approximated by a momentum-independent, but k- 
dependent running gauge coupling. A similar truncation is used for the ghost con- 
tribution. In sect. 5 we truncate, in addition, the most general form of the four- 
and six-fermion interactions. 

We are interested in the evolution of the two- and four-point functions for the 
quarks. The respective flow equations for these quantities obtain by taking the 
second and fourth functional derivative of eq. ( |2.33| ) at ip = ijj = 0. We label the 
different contributions on the r.h.s. of eq. (|2.33|) by 



is the standard contribution of a pure fermionic theory. The remaining terms 7^^, 7^ 
and 7c involve R^/f^ and reflect the contributions from gluons, whereas e gives the 
ghost contribution. The term 




(2.34) 



and discuss them separately. The first term 




(2.35) 




(2.36) 
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involves only a trace over gluonic degrees of freedom and accounts for the contribu- 
tion of gluon fluctuations around the ■^/'-dependent classical solution. The contribu- 
tion of 7a to the fermionic two- and four-point functions is given by the dependence 
of r[,^''['?/', 0]"^/ on ip. Relevant contributions to 7^ therefore arise from terms 

in Tk[ip, A, 0] which are either quadratic in A and also depend on or are cubic or 
higher-order in A. 

We use a truncation where the first sort of terms is absent and no relevant con- 
tribution to 7a would be present for an abelian gauge theory. For nonabelian gauge 
theories we get contributions /^from the three- and four-gluon vertices in [■?/', A]. 
We approximate here these vertices by the (standard) lowest order expressions which 
are obtained from functional derivatives of F^^F^'^ . More precisely, we use on the 
r.h.s. of the flow equation 

^ (rfl^, = 01):; (.,.') = {v\a\):j(^-^') (2.37) 

with 

Z-/ {V[A])1'^ = -{D'[A])yjl^ + 2i~g {T^f ^ + {bM]D^[A])\ (2.38) 

Here -D^[A] = — igA^^Tz represents the covariant derivative in the adjoint repre- 
sentation with gauge coupling g and F^^^, is the nonabelian field strength associated 
to the gauge field A^. In this truncation the three- and four-gluon vertices are 
parametrized by two running parameters g{k) and Zp{k). The running renormal- 
ized gauge coupling g^ is related to them by 

gl = fz-/ (2.39) 

We observe that A, 0] is invariant under global gauge transformations of ip and 
A. The expression for 7^ does not explicitly depend on ip in our truncation and 
7a[^ = 0, y4fe] or rfc[0,y4fc,0] cannot contain a term linear in Ak. There is therefore 
no contribution from 7^ to the flow equation of the fermionic two-point function. An 
estimate of the contribution to the four-quark interaction from 7^ therefore amounts 
to a computation of the gluon contribution to the evolution of the term quadratic 
in A in Tk[ip = 0, A, A = 0]. 

Similarly the ghost contribution e is (with the approximation (|2.3CI|) ) only a 
functional of A, containing terms quadratic in A (and higher orders). We also 
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observe that 7a and e only account for the gluon and ghost contributions to the 
effective gluon propagator, whereas the contribution from quark loops is implicitly 
contained in 7^. We note that the latter is not distinguished any more from any other 
fermionic contributions, as, for example, from an explicit four-quark interaction in 
Ck[^]. The contribution 

7c = 2^L' (^^j At (2.40) 

describes the effect of the "classical" change in the infrared cutoff as k is lowered. 
It is quadratic in the classical solution and therefore gives a contribution to 

the fermionic four-point function. Finally the piece 7^,/, involves the explicit ip- 
dependence of the classical solution ~ dAj^/dip. It contributes to the running of the 
fermionic two- and four-point functions. 



3 Heavy quark approximation 

In the limit of infinitely large quark masses our formalism simplifies considerably. For 
euclidean momenta we can omit in eq. ( p. 33 ) the terms involving the inverse fermion 



propag ator (rl^^[7/;]+i??Vi since their contribution is suppressed by inverse powers 
of the quark masses. In the language of the last section this results in 7^ = 0, 7^,/, = 
0. The remaining computation amounts to an investigation of the pure gluon theory 
with static quarks. This is done most easily in the language where the gluon fields 
are kept explicitly and the relevant effective action is T^lip, A,A = 0]. If one wants 
to extract the effective four-quark interaction, one needs the /c-dependent effective 
gluon propagator and the effective vertex tptpA. We first describe for arbitrary quark 
masses the general framework how an effective four-quark interaction obtains from 
"gluon exchange" in the formulation where both quark and gluon degrees of freedom 
are kept explicitly. We then specialize to the heavy quark limit. 
Let us consider in [■?/', A, A] the term quadratic in A 

= IJ (0i^;(-9)(rr[^ = 0, A = 0, A = 0])T^A^,{q) (3.1) 
and parametrize the most general inverse gluon propagator by 

;Ff [ij = o,A = o,A = 0])'' (q) = {GA{qn + ^^A(g)g.g'^)5f (3.2) 
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For the quark-gluon vertex 

J- I 



= J J^.^Mp)G4p^ qmT.Vm + q)Al{-q) (3.3) 

we make the approximation that is a simple function not involving Dirac matri- 
ces. Knowledge of G^, Ha and permits to compute the classical solution in 
order ipip 

,(0). r.u,_^ f d^P 



{AT{q)): = -S;iq) J j^,G^{pAWa{pmT^\'r,{p + q) (3.4) 
Here S = (T^^O] + Rk)'^ is the gluon propagator in presence of the infrared cutoff 

(Rl^'^y/Aq) = (RkiqK + Rk{q)quqn8l (3.5) 

and reads 

S;{q) = (GAiq) + Rk{q)r' {^^ - q'^q.iHAiq) + Mq))' 

[GAiq) + Rk{q) + q\HA{q) + Rk{q))V} (3-6) 

Inserting the classical solution into (|3^ ) and ( p.3|) and accounting for the term A^'^^S' 
(2.22) we find the effective quark four point function 

(^) 1 r d^p d^p' d^q 



(2^ (2^ (2^^-'^^)^^^^' 



mp)r{T;),^rAp + q)}mp')iuiT^)M{p' - q)} 



1 [ d^pi d'^pi 4 

7 72^-72^^'^^ %1+P2-P3-P4) 



with 



and 



2 J {2TxY"'{2Tiy 

{FlipuP2,P3,P4)M{pi,P2,P3,P4) + ^2(Pl,P2,P3,P4)A/'(pi,P2,P3,P4)} (3.7) 



-^1 = G^{-pi,Pi - P3)G^{p4,,P2 - P4)iGAiPl - Ps) + Rk{Pl - P3))~'^ 

F2 = G^{-pi,pi~p^)G^{pi,p2-pA){HA{pi-pz) + Rk{Pi-Pz)) 
{.Ga{,Pi -ps) + Rki.Pi - P-i)y^[GA{.Pi -Ps) + Rki.Pi -P3) 

+(pi -p3fiHAipi-p3)+RkiPi-Ps)T' (3.8) 

^fiPuP2,P3,P4) = m-Pl){A-h)iT^),'i^]i-p3)} 

mP4){i>2-k)iT^)M{P2)} (3.9) 
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The curled brackets indicate contractions over not explicitly written indices (here 
spinor indices), k,i = l...Nc are the colour indices and a,b = l-.-Nf the flavour 
indices of the quarks. By an appropriate Fierz transformation and using the identity 

iTl/{T.),' = \S',61 - ^^6i6i (3.11) 

we can split into three terms [0] 

M = M^ + Mp + Mp (3.12) 
= -^{V^:(-Pi)^.^(P2)}{^^(P4)^^J(-P3)} 

l {^:i-Pih'i^1iP2)} {4iP4h'i^]i-P3)} (3.13) 

I {^:i-pih,i^1{p2} {i^l{P4h'^]i-P3)} 

\ [i'i{-Pi)i,i'^'{P2)] {i^'Mi'i^r^i-P^)] (3.14) 

Mp = -^{V5:(-Pi)7mC(-P3)}{V'^(P4)7'^^?(P2)} (3.15) 
In terms of the Lorentz invariants 

s = (p^ +p2)^ = (P3 +P4)^ 

t = {Pl-P3f = {p2-P4)' (3.16) 

we recognize that the quantum numbers of the fermion bilinears in A^^r correspond 
to colour singlet, flavour non-singlet scalars in the s-channel and similarly for spin- 
one mesons for Aip. In analogy to ref. we associate these terms with the scalar 
mesons of the linear a-model and with the p-mesons. The bilinears in the last 
term Aip correspond to a colour and flavour singlet spin-one boson in the t-channel. 
These are the quantum numbers of the pomeron. We observe that in the heavy 
quark approximation where T^^l arises only from "gluon exchange" , the coefficients 
of the quark interactions in the a, p and pomeron channel (|3.12| ) are all given by the 
same function Fi. 

The general quark bilinear is conveniently parametrized by the real functions 
Z^{q) and rha{q) 

= E / ^M<l){ZMYq, + ma{qmm) (3.17) 
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The fc-dependence of the functions Ga, Ha, G^, and rha relevant for the two- and 
four-point functions for the quarks can now be studied using the evolution equation 
( p.28| ) for rk[ip,A,A = 0]. In the truncation where only the terms (|3.2|) , (|3.3| ) and 
( p.lTp are kept, it is easy to see that the contributons to the ^-dependence of G^, 
and rha all involve quark propagators. In the heavy quark limit they can therefore 
be neglected for euclidean external momenta. Only the fc-dependence of r\f'2 needs 
to be considered. For Z^ and G^ we may take appropriate momentum-independent 
"short-distance couplings" 

Zi,{q) = 1 
~ i 

GAp^Q) = Z^{m^)g{m^) = g{m^) (3.18) 

with renormalized gauge coupling g taken at the scale k = and the heavy 
quark mass. We also may identify k = with the "ultraviolet cutoff" or the scale 
where the initial values for the flow equation are specified, i.e. 

Zpim^) = 1 
GA{q]k = m^) = (3.19) 

Solving the flow equation for GA{q) for — ^ yields the effective four-quark in- 
teractions for momenta much smaller than the quark mass. For = (see next 
section) the effective four-quark interaction is fully determined by 

To? = —g^im^) J -^^--^^W^^iPl + P2 - P3 - Pi) 

\imG^^{pi -Pa) ^M{puP2,P3,P4) + j^^^^^^^{PuP2,P3,P4)^ (3.20) 

We finally should mention that the heavy quark potential or the scattering am- 
plitude for heavy quarks cannot be extracted directly from the four-point function 
at small momenta <^ m^. For these purposes the momenta appearing in 
should be taken on-shell, i.e. pf = P2 = P3 = P4 = Their size is therefore 

not small as compared to m^. For on-shell momenta the vertex function G^{p,q) 
becomes a function of q^. The flow of G^{q) does not vanish for k"^ ^ and 
we need to supplement the computation of the /c-dependence of Ga(q') by a corre- 
sponding one of G^{q). The heavy quark pontential can then be extracted as the 
(three-dimensional) Fourier transform of F{q) = G'^{q)G^^{q). 
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4 Scale dependence of the gluon propagator 



In this section we compute a flow equation for tlie scale dependence of the effective 
gluon propagator in the pure Yang-Mills theory. As explained in the last section, this 
is the central piece which determines the heavy quark interactions. We will derive the 
flow equation for Ga in the approximation that the vertices can be extracted from a 
term ^ZpF^^F^^" in the average action, with running wave function renormalization 
Zp and running gauge coupling. The resulting functional form ([4.9|) of the r.h.s. 
of the flow equation corresponds to 7^ — e and is valid beyond the heavy quark 
approximation. The second part of this section discusses qualitative properties of 
the solution of the flow equations. This part is valid only for the pure Yang-Mills 
theory or for the heavy quark theory. For the derivation of the evolution equation 
we keep the most general gluon propagator. The only truncation concerns the 
momentum dependence of the effective three-gluon and four-gluon vertices and the 
ghost sector. One finds that the /c-dependence of the functions Ga{(i) and HA^q) is 
governed by the flow equations 



d_ 
di 



GA{q) = N^glZp 



dtl{GA{q')+Rk{q')y Zp 



1 



-- (GA(g') + {GA{q + g') + Rk{q 



q')) 



-1 



72 



5g2 + 2 
2 



^ 16 ,2 



10 {qq 



A 2 



biq') [ 2/ + 10 



1, 



q 



q 

/\2' 



.f\3 



+ lp,^\q')P,^\q + q') 



q 

{qq'? 
q2 



i2 



/\2' 



(4.1) 



and 



-^N^dZp 



(Ta' ~ 
aq 



(2vr) 



{GA{q') + Rk{q')y' {GA{q + q') + Rk{q + g'))"' Z] 



"40 {qq 



'\2 



2\2 



10 



,/2 



+ 10 



{qq') 
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U{qq'f Iq'^ 4: {qq'f 8 {qq 
y g2g'2 + 3^ " 3 (g2)2 



/\3 



3 (g2)2g'2 



3 (g + g')^ 



■-P-hi(g')Pgh'(? + 9') 



\qq 



A 2 



V2 



+ 4 



r2l 



Here we use 



(g2)2 



(4.2) 



(4.3) 



GAiq) + Rkiq) + (i^A(g) + RkiqW 
and note that the partial derivative dt acts only on the explicit infrared cutoff terms 
Rk and Rk- The parts involving the effective ghost propagator arise from the 
ghost contribution ~ — e. After performing the q' integration the evolution equations 



?T]), (^4.2| ) can be interpreted as two coupled nonlinear partial differential equations 
for the functions Ga,Ha which depend on two variables k and g^. They describe 
the scale dependence of the gluon propagator in the approximation where both 
the three-gluon and the four-gluon vertex are given by a single renormalized gauge 
coupling Qk and similarly for the ghost gluon vertex. 

Even in this approximation the flow equations are lengthy and difficult to solve. 
A simplification occurs if we take for the gauge-fixing term 



1 

2a 



d^x{d^Al 



(4.4) 



the gauge parameter a — > 0. In this limit Ha diverges ~ ^ and b{q) approaches one. 
We may define a /c-dependent renormalized gauge fixing parameter an by0 



HAiO) 



-1]Z 



(4.5) 



with limfc-^oo Zp = I and limfc^oo or = oi. The evolution equation for aR follows 
from 



d 



(4.6) 

0. We conclude that 



dt ^ " dt 
where we note that ^HA{q) has a well defined limit for q^ 
aR = is a fixpoint which is infrared stable for din Zp/dt > 0. The approximation 
a/j — s> remains therefore stable in the course of the evolution. We further observe 

•^We choose Rkq^ = — 1^ Rk such that b{q) = 1 + 0(q;_r) for all values of q'^. 
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that for aR ^ the function F2 defined in eq. ( p. 81 ) equals —ps)^ and only the 

function GA{q) determines the effective four-quark interaction. The flow equation 
for Ga can be written in a more compact form using 



PA{q) = Z]^\GA{q) + Rk{q)) 
We also approximate the ghost part of T^. ^ by 



PM=PA{q), dtPM=dtPA{q) 



(4.7) 



(4.^ 



This yields for an ^ 



d 



j-GAiq) = N^glZp 



6 



13g2 - 14(gg') + lOq'^ 



-10 



{qq 



l\2 
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{qq 



A2 



V2 



-4 



{qq' 



/\3 



^ q"- q^q'"^ ' g'^ (g + q'Y 

The evolution equation ( [4. 91 ) is the central equation of this section. It is a partial 
nonlinear differential equation for G^(g^; k) which can be solved numerically. We 
observe that the r.h.s. of the flow equation ( |4.9|) involves not only Ga('?) but also 
the renormalized gauge coupling and the gluon wave function renormalization 
constant Zp. The precise definition of these /c-dependent constants is a somewhat 
subtle issue. We could define g and Zp (and in consequence gk) in terms of the 
effective three-gluon vertex ~ gZpA^dA and four-gluon vertex ~ g'^ZpA'^ which 
enter the approximation for T\. used in (2.38). In this way g and Zp are expressed 
in terms of third and fourth functional derivatives of Tklip = 0, A, A = 0] evalu- 
ated at A = and projected on the appropriate index structures. One also has to 
choose appropriate momenta for the external legs for the effective vertices, as for 
example the limit where all momenta approach zero. Evolution equations for the 
scale dependence of Zp,g and gk could then be computed from appropriate func- 
tional derivatives of the flow equation (|2.28|). We will use here a simplification and 



+ 



q q q' — {qq 
v2' 



l\2 



(4.9) 



GA{q)=rni + G^^q' + G^hq'f + ... 

and identify Zp with G^^ . There is an obvious limitation to this approximation since 
G^'' may turn negative for small k and any reasonable choice of a wave function 



approximate Zp by the coefficient of the q^ term in Ga- More precisely, we expand 
for small q^ 

(4.10) 
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renormalization requires positive Zp. For the flow equation for tlie renormalized 
gauge coupling g^^ we rely on the fact that the first two coefficients of the /3 function 

^""^ - = A - c,jj\^ - ... (4.11) 



are universal 



dt 167r2 "CieTT 



22N, 204 ^ . 



ci = ^ C2 = —K (4.12) 

In the region of large gt we may also use nonperturbative estimates of (3g2 derived 
by different methods 0. An ansatz for (3g2 combined with an estimate of fjF fixes 
also the evolution of and provides all information needed for a numerical solution 
of the flow equation. 

We concentrate first on an analytic discussion of a few prominent features of the 
solution of equation ( [4.9| ). The evolution equations for the mass term fh\ and for 
G^'' are easily derived by expanding the r.h.s. of eq. ( |4.9| ) in powers of q. One 
obtains 

|mi = N^gfZp J 0^^d, [lPX\q') - ^g'^P/(g')} (4-13) 
and the renormalized dimensionless mass term 

m\ = m\Zp^k-^ (4.14) 

therefore obeys 



^fi\ = (-2 + f^p)mi - ^gt (^-/^i - fX2 ) (4.15) 



V2 



Here we have defined the integrals, with x = q' 

1 roo , 

^ln = -ofc"/ dxx-^-^dtPX^{x) (4.16) 

2 JO 



The evolution of the ration m\/ g^^ 



2 



is characterized for small g1 by an approximate infrared unstable fixpoint 

gl 1287r2 



(4.18) 
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Indeed, we can use for small g1 the lowest order expressions 

GA{q) = Zpq^ 
dtPAiq) = ^P(g) (4.19) 

with 

P{q) = q' + Z-/R,{q) = ^-^^-^ (4.20) 

l-exp^-f^j 

such that l\ i = ^1 = ^Ja,2 = ^2 = |- We also neglect Pg/g^ — fjp as compared 
to two. For small gf. the general identities for the dependence of [■?/', A, A] on 
the background field A 0, or, similarly, the generalized Slavnov- Taylor identities 
imply that the fc-dependent mass term is indeed described by this fixpoint 
(cf. Appendix). We conclude that for small gt the mass term induces only a small 
correction in the momentum-independent part of P4 



P + Zp^m\ + 0{q'') 

kH^-^.9i]+0{q^) (4.21) 



3iVc ,^2 
1287r2^^ , 

The behaviour of ffi\ near the confinement scale where g1 becomes large is more 
complicated and best described by evaluating directly the relevant identities It 
is conceivable that a negative r.h.s. of ( [4 .131) drives ffi\ to a positive value for A; ^ 0, 
but we find this scenario not very likely. 

The flow equation for G^^ can be written in the form 



^^(1) _ J^.27_^Qi;4 p,_6 

dt^ 967r2 



G7 = ^giZFi3llX2-5m'A 

y'^^^^ievr 



NcCAT^ (4.22) 



where we have defined the integrals 

<n = -^k'-'-' I dxxW, {Pl{x)Pr{x)) (4.23) 
with Pa = -^Pa- We note that the r.h.s. of eq. (f4.22|) is positive for positive ca 



CA = ^(31/1,2 -5m^,4) 



If <<..a.{p-(l-A.'P^P-)} (4.24) 
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Let us first consider small values of the gauge coupling where we can approximate 
BiPa = ^P- In this limit we obtain ca = 10- For the definition 

Zp = G^}^ (4.25) 

and therefore finds 



riF = -^InZp = -—N,CA^ (4.26) 



^I^Zf = --N,ca-^, 
dt 3 167r2 

For a solution of (^4.26 ) it is convenient to compare fjp with 

where 6a may depend on or k and reflects the deviation ibom. the one-loop (3- 
function for which 6^ = 1- As long as the fc-dependence of the ratio hA/cA can be 
neglected, one obtains 

Mk) (g\k,)V UCA 

ZAko) \9'ik)) ' ^ 226a ■ ^ ' ' 

In the one-loop approximation decreases logarithmically and reaches zero at 
the confinement scale Aconf- We conclude that this is exactly the scale where Zp 
would vanish. In this language the divergence of the renormalized gauge coupling is 
actually entirely due to the vanishing of Zp: From (|4.28|) one obtains 

~9\k) ^ (g\k)y-' 
~9\ko) Wiko)) 
ZF{k)g\k) fg\k)Y-'' 



Zpikoyg^iko) Wiko) 
Zp{k)g{k) ( g{k) ' 



1-37 

(4.29) 



Zp{ko)g{ko) \9{ko)J 
and for 7 > | (cf. (|4.28| )) both the unrenormalized three-point vertex ~ Zpg and 
four-point vertex ~ Zpg^ vanish at the confinement scale. 

As mentioned before we should not use ( [4.25| ) for k in the vicinity of the con- 
finement scale. We therefore propose [|1^ to keep Zp independent of k for k < k^p 
where knp is defined by 

2 47r^ 

^We have also computed ^Ha(S^) — — 4^5^ (8^^ 2 + S'^A 4) ^'^'^ ^'^'^ ^^''^ small 5^ that 
^ {^Ha{Q) + G^'^ vanishes as required by the Slavnov- Taylor identity for perturbative trans- 
versality. 
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The definition 

Zf= { ^"^^ ' 4.31 
\ Zf{K,) for gl>'^ 

allows to separate the issue of vanishing from the choice of the infrared cutoff 
(i.e.Zi;'). 

The vanishing of the term G^'*g^ would have important consequences for the 
behaviour of the propagator ~ G^^. We should therefore investigate if this feature 
is likely to survive beyond the approximation of small gauge coupling. Turning back 
to eq. ( [4. 221) we observe that G^"* can only remain positive for k —>■ if turns to 



zero or if ca vanishes or becomes negative. In view of eq. (|4.29|) the first alternative 



seems not very likely. (This probably generalizes if we go beyond the approximation 
leading to ( [4.1| ) and consider general momentum- dependent three- and four-point 



functions. Then in eq. (f4.22|) g has to be replaced by an appropriate momentum- 



weighted average of these vertices.) In order to investigate the sign of ca we write 
the integral ( |4.24| ) explicitly as 



c. = H ,,,p-^s + (4.32) 



with 



26 Jo ^ 31 VP. / 31 PaS 



S{x)=dtPA{x) = ^^P{x)-f,F{P{x)~x) 



2^ - VF ] {P{x) - x) (4.33) 



The integral is dominated by the region x ~ /c^ and can turn negative only if the 
bracket is negative in this region. As an illustration we take Pa = P + itl\ + 
Kx^, m\ = Zp^Th\. This implies 

xPa Px + 2kx'^ 



Pa P + rriA + nx"^ 



(4.34) 



to be compared with 



and 



xS _ xP / fiFk^\ ^ P 



1 - -^f- 1 + 1 - (4-36) 
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We conclude that negative m\ (cf. ( [4.18| )) and positive k tend to lower ca- A 
vanishing ca for large k, cannot be excluded without a more detailed investigation. 
On the other side, if the integral relation 

—g\k)cA{k) = — — 4.37 

Aconf k 13iVc 

could be fulfilled for Aco„/ > 0, the coefficient would vanish at the confine- 
ment scale Aconf and presumably becomes negative for k < Aconf- We emphasize 
that a negative value of Th\ or of G^'' for small k would independently indicate 
that the groundstate does not correspond to the perturbative ground state ^4^ = 



Even before reaching the confinement scale, the gluon propagator has to be 
modified: Whenever the term G^''g^ becomes comparable to G^^\q'^y for ^ /c^, 
i.e. for ~ G^'^k'^, the gluon propagator cannot be approximated any more by 
the inverse of g^! As an example for a plausible form one may consider 

GA{q) -m\ + G^^g^ + ZpK^^^ (4.38) 

1 + og^ 

where 

' ^"'^ (4.39) 



is determined by the requirement that for large g^ one expects GA(g) = g^ indepen- 
dent of k. (This holds up to neglected logarithmic corrections.) Assuming that for 
k ~ knp the term m^ + G^^g^ is small as compared to the {q^Y term the approximate 
form of the propagator 

would be close to the one corresponding to a confining potential. 

It is obviously difficult to find an analytical answer to all these questions and 
it seems preferable to solve the fiow equation ( [4.9| ) numerically. A numerical inves- 
tigation has been performed by B. Bergerhoff and the author [Q. We show here 
only a few first results. In fig. 1 we plot the scale dependence of the wave function 
renormalization Zp as defined by 

d 
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with Aqcd the two- loop confinement scale. This definition implements the idea that 
Zp should be taken constant for very small k (cf. eq. (4.31)) in a smooth way. One 
observes the decrease of Zp as k is lowered corresponding to eq. ( [4.28|) . We have 
started the running in the perturbative region at A; = 40 GeV with G'A(Q')|fc=4o Gev 
given by the one-loop perturbative expression (containing the infrared cutoff). The 
latter was normalized by SGa/ dq^\q2^Q = 1, which explains the starting value of Zp 
(40 GeV) somewhat larger than one for the definition ( [4.41| ). In order to concentrate 
on the deviation of Ga(q') from the linear dependence on we introduce the quantity 

Here the mass term G'a(O) is subtracted from the inverse propagator and the lead- 
ing perturbative g^-dependence is divided out. In the classical approximation the 
expression in the bracket equals one, and a nonvanishing value of xIq) is entirely 
due to quantum fiuctuations. Since a computation of x{q) involves a numerical 
derivative of a small difference, one needs a numerical solution of the flow equation 
for Ga{(i) with a relatively high precision. We observe that %(<?) plays the role of 
a momentum-dependent anomalous dimension. Within renormalization-group im- 
proved perturbation theory one expects for /c = (compare eq. (4.26)) 

XM = f iV.f^ (4.43) 

where (7^(/i) is the running gauge coupling at the scale fi. In flgs. 2 and 3 we compare 
the numerical determination of xIq) with the renormalization-group improved one- 
loop perturbative result (4.44), with g the two-loop running gauge coupling. Up to 
a scaling factor of about 10 % the two curves asymptotically coincide for large as 
k goes to zero. This is not a trivial result since no assumption of this type enters 
the flow equation ( |4.9| ). We note that a propagator ~ for small corresponds 
to limfc^olinig2^oX(Q') = 1- ^ a behaviour extends effectively over a certain 
momentum range xIq) should develop a plateau at one in this range. As k is lowered 
we see in flg. 3 a sizeable increase of xIq) for small q^. We also observe a tendency 
of an extension and flattening of the maximum somewhat below one. This tendency 
should stabilize as k goes to zero. We have stopped the running at A; = 400 MeV 
since for small k"^ and g^ the approximations leading to (4.9) become doubtful. At 
least part of the momentum dependence of the gluon vertices according to (A. 55) 
should presumably be included. 
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5 Flow equations for light quarks 



The formalism for integrating out the gluon degrees of freedom needs not to be 
restricted to the heavy quark approximation. We want to derive in this section the 
general flow equation for the quark two-point and four-point function corresponding 
to (2.33). For simplicity we mainly consider Np massless quarks - the inclusion of 
mass terms is straightforward - and we work in the gauge with a = 0. At some 
appropriate short distance scale we start with the "classical action" (cf. (3.7)) 

Tkli-] = Z^,j - Izjsl . (5,1) 

The flow equation then permits to study how changes its form as k is lowered. In 
particular one is interested in pole-like structures in the quark four-point function 
which would indicate the formation of meson bound states [p!5[| ,[3. 

In order to establish the flow equation we have to collect various pieces which 
have been discussed in the previous sections. We begin with the contribution from 
the gluon and ghost fluctuations 7yi ~ ^ + 7c which only contribute to the four-quark 
interaction Ffc^4 [■?/;]. The direct contributions from gluon and ghost loops read 

7A - e = ^ / ^ [Af{q))l {Af{-q))[\^^GA{qK + |^A(g)g%} (5.2) 

The functions ^Ga{(1) and -^Ha{q) have been computed in the last section (cf. (|4.1|), 
( [4.2|) , or (^^) ), where we have indicated by a hat that only gluon and ghost contrib- 
utons should be included here, in analogy to the heavy quark approximation. The 
classical solution A^^^ is given by eqs. ( p.4|) and ( p.6|) , where Ga and Ha characterize 
now the gluon propagator without reference to the heavy quark approximation. We 
consider again the limit a — * where S^{q) = {G a + Rkl'^i^^^ — q^qiil q^)- Including 
also the contribution 7^ which is similar in structure we obtain 

7A-e + 7c = -2j/ (^•••(^(2v^)5(p,+p,-p3-P4)• 
/l(Pl,P2,P3,P4)(A^+^^^^4^^) (^-2) 

with 

fiiPi,P2,P3,P4) = -G^{-pi,pi -p3)G^{p4,p2 - Pa) ■ 
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/ 9 - d \ 

{GA{pi-Pz)+Rk{pi-Pz)y^ [-q-Ga{Pi-P^) + qiMpi -Ps) 1 (5.4) 

(In the heavy quark approximation there is no difference between Ga and Ga- Also 



G^ is independent of k such that /i = (|3.8| ).) The function /i involves the quark 
gluon vertex G^ and the gluon propagator Ga which also enters in the determination 
of §-fiA. 

Up to this point the only approximations involve the three- and four-gluon ver- 
tices entering ^Ga and the ghost sector. They have been discussed in sect. 4. We 
may in addition also truncate the quark gluon vertex and use for light quarks the 
ansatz 

G^{p,q) = Z]l'^Z^gk (5.5) 



The running of the renormalized gauge coupling gk is now determined by the (3- 
function including quark contributions and Zp may be identified with G^^ for k > 
knp (cf. sect. 4). the lowest order truncation for Ga{(i) would be 



GA{q) = Zpq' (5.6) 

such that -Pa(q') is replaced by P{q) in the equation ( |4.9| ) for ^Ga- The detailed 
discussion of the last section shows, however, that this approximation becomes in- 
valid for k in the vicinity of the confinement scale. There one should rather use 
a truncation of the form (|4.38| ) or similar. The corresponding flow equations for 
Th\,G^^^ and k include now additional contributions from quark fluctuations and 
have to be computed from the evolution equation for ^GA{q) in the formulation 
where both quark and gluon degrees of freedom are present. Fortunately, the for- 
mation of meson-bound states occurs typically at a scale substantially higher than 
the confinement scale. This gives the hope that important features of meson physics 
can already be extracted using the truncation ( |5.6|) on the r.h.s. of the flow equation 
and do not need a very detailed understanding of gluon condensation phenomena. 

With the truncation ( |5.5| ), ( |5.6| ) the function /i ( |5.4| ) only depends on the Man- 
delstam variable t = {pi — Ps)"^ 



1 dt 



^P{pi - ps) - fiF {P{pi - P3) - {Pi - Ps)') + N.glGipi - p3) } (5.7) 
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where 



^G^(g) = N^glZpGiq) (5i 



is given by the r.h.s. of ( |4.9|) with Pa{(1) = and dtPAio) = ^P{(l) ~Vf{P{(i) — 
q^). The first two contributions are proportional dtP~^ and simply reflect the change 
of the infrared cutoff contained in P~^ in the "classical action " ( p. 11 ). Only the last 
term ~ Q describes how additional quantum fluctuations are included as k is lowered 
- in this case the gluon and ghost contributions to the gluon propagator. If we omit 
the contribution ~ fip in dtPA^q) the funtion Q{q) reads explicitly 



13,^ - 14(to') + 10,'= - 10^^ - 22^£ 



q" q' 



{qq'f q^ q^q'^ - {qq'f 



(5.9) 



We observe that even for our simple truncations the function G{q) has a compli- 
cated momentum dependence. A solution of the flow equation for the four-quark 
interaction will go far beyond the effects of a running gauge coupling in the one-loop 
approximation . 

The term contributes to the flow equation for the two-point and the four- 
point function. The contribution to the two-point function can be extracted from 
(2.33) using the classical solution (p.4|). With a = and G^{p,q) = Z^Zp'^g^ the 
lowest order classical solution reads 

{Af\q)r, = -{GA{q)+Rk{q))-'Z^Zl/'g, 

'^'^-Mp) (7^- ^1 [T-^VrAp + q] (5.10) 



(27r)4""'"'V' q 
and the next to leading contribution A^jl^ obtains as 

[A^\<l))[ = -^9lZ-,"'zlfrPA\<D 

^(Pl - - P2 + P4) ({V5(-Pl)7i.7;^(-P3)}{^(P4)7''^«;V'(P2)} 
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{Pl-P3)u{P2-PiT 
{Pl -P3)^(P2 -Pif 



{i^{-PlWl -^z)Tyi^{-P^)}{'^{PiW2 -A)Twi^{P2)} 



{P2-P4) 



;{^{-Pl)(]^2 -A)Tyij{-P3)}{^{P^)¥2 -A)Tw^{P2)} 



{^(P4)7'^T^^(P2)} 



(5.11) 



With the ansatz of a flavour diagonal kinetic term and mass term 

r(')[0] = Z^(c„(gy + m„(g)7')5,»5i(27r)^5(g - 4) (5-12) 
one finds the following contribution to the fiow equation for the two-point function 



(2) _ K 



2N, 



a 

[(ca(g)+rf (9))(^+2^^^(^+|i))-3m„(g)f]^„(p)} (5.13) 

The contribution to the four-point function 7^^ can be obtained similarly using 
(5.11). 

For a computation of the purely fermionic contribution 7^ we will use the fol- 
lowing truncation for the term quartic in the fermionic fields 

{K{Pl,P2,P-i,PA)Ma + \p{pUp2,Pi,PA)Mp 

+>^p{Pl:P2:P3:P4)Mp + Ki.Pl, P2, Pz, Pa)M} (5.14) 

This yields a contribution to the fiow equation for the two-point function 

2 

2iVcAp(-g, q, -p,p)4 - —Xp{-q, q, -p,p)4 

jy2 _ 

' -K(-q,q,-P,p) ((p^-q^)<i+'^(q^-{pq) 
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[ca{q) + rf ^(g)) gV [8iVcA.(-g, g, -p,p) 

-J-Ap(-g,g, - ^'^^^ ^\ q-pf>^n{-q,q, 



^ (c.(g)+rf (g))V-»^^(g) [2A.(-g,p,-g,p) 
+2Ap(-g,p, - 4Ap(-g,p, . 



(5.15) 



For j^^-^ we present here only the case Ap = Ap = An = in the chiral hmit (m^ = 0). 
Omitting all contributions except "jlp^ one finds 

d 

g^K{Pl,P2,P3,P4) = '2Vi'K{Pl,P2,P3,Pi) 



d^q qM-Pi-p-i) 



Hq) + rt\q)]-' 



^^^^ ' {2-kY qi[q-p^-p^Y 
dt[c{q-pi -P2) +r'j^\q-pi -^2)]"^ 

Kipi,P2,q, -q + Pi + P2)Kiq, -q + Pi + P2,P3,P4) 



(5.16) 



For the special case (c(g) + rl^\q))^^ = exp — exp (— fj) ,Vip = ^ this re- 

produces the fiow equation of ref. 0. The r.h.s. of this equation should now be 
supplemented by the contributions from 7a — e + 7c + 7^^/; which have been discussed 
before. Furthermore, a better approximation should include the contributions from 
Ap etc. 

Besides this, eq. ( |5.16 ) involves the explicit momentum dependence of the 
fermion kinetic term, i.e. the function c(g). The scale dependence of the quark 
propagator can be computed from 7^^"* and 7^^. Combining eq. ( ^.15|) with ( ^.13|) 
and (|5.12|) yields the fiow equation for the kinetic term 



d 



Z^Ca{p)) = 



dt 



d^q (N^-l 



(2vr) 



9k i^rkip - g) - VFTkip - g) 



dt 



{cM + r'if\q)y q' + mliq) 



(c.(g)+rf(g)) (2g^ - 3(,g) - + 4M 



^rf^(g) -r/^rf^(g) 



(ca(g)+rf (g))'g^ + m^(g) 



1 -2 
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{ca{q) + r'^\q)y q"^ - ml{q) 



2 \ (pq) 

2NcXp{-q, q, -p,p) - —Xp{-q, q, -p,p) ) 



p^ 



-K{-q,q,-P,P) 2g - [pq) 



(pq) ( d 



p. 



^^(c,(g)+rf (g))V-m,2(g) 



[2Xa{-q,p, -q,p) + 2Xp{-q,p, -q,p) - 4Ap(-g,p, -q,p)] \ 



(5.17) 



The right-hand side of this equation involves the gauge couphng and the effec- 
tive inverse gauge field propagator Pj^ as well as the effective fermionic four-point 
vertices A^, Ap, Xp and A„. It is instructive to study this equation in the "classical 
approximation" for the four-quark vertices A^, i.e. 



K{Pl,P2,P3,Pi) = Ap(pi,P2,P3,P4) = Ap(pi,P2,P3,P4) 
= Id {PX\Pl - Ps) + PX\P2 - P4)) 



(5.18) 



K{Pl,P2,P3,P4) 

Using 



:9l 



+ 



, (Pi - P3)^Pa{Pi - Pz) {P2 - P4)'^Pa{P2 - P4) , 



dtPX\q) = -q'PX\q)[^^rk(q)-f,Fn(q)^ 

{ca{q) + r'j^\q)fq^ - ml{q) 



Ca{q) + rf\q)Yq^ + rnl{q) 



' d_ 

at 



(5.19) 



one obtains 

d 



dt 



Zi,Ca{p)) 



oAr 9k 



2Nr 



A /2(pg) 2yV-(w)(/ + q2)' 



(27r)^ 



p^{p - qf 



dt{PX\p-qy 



Ca{q)+rf\q) 



(c.(?)+rf (?))V + m^(g), 



(5.20) 



We observe that this expression corresponds to the formal dt derivative of the stan- 
dard one-loop correction to the fermion kinetic term in presence of an infrared cutoff 
in the propagator. 
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The fermionic wave function renormalization can be defined by Ca(po) = 1 
for suitable po and a. We will choose here po = and use Ca corresponding to a light 
quark. Defining the anomalous dimension 

d 

Vi^ = -Q^^^^i^ (5-21) 

one obtains for rria = 
3 - 1 " 



4 iV. "^^ / ■^.^29t{ic{q)+r^^\q)r'PX\q){l-q'pA{q)PA\q))} 

(5.22) 

where Pa = dPA/dq^. In the perturbative limit c(g) = 1, PA^q) = P{.q), dtPA{,q) = 
■^P{q) this yields a vanishing fermionic dimension 

r/v- = (5.23) 



We also may extract the anomalous mass dimension by expanding eqs. ( p. 13 ) 
and ( |5.17| ) in linear order in rria 

d 

— m«(0) = UmrriaiO) (5.24) 



One obtains 



■J{2nrm{0)\2 ' \c{q) + rf {q)W' ^ ^ 

+^'77\ 1), , [4iVcA.(-g, q, 0, 0) - ^A,(-g, g, 0, 0) 

iV^ - 1 



g^A„(-g,g,0,0)] 



which reduces in the "classical approximation" Ao-(— g, 0, 0) = Ap(— g, g,0,0) 
9feAn(-g, g,0,0) = \g^PX^{q) and for ma(g) = ma(0) to 



= gi / -^-Tjdt TT^ (5.26) 

2 AT, ^'M2vr)4 *(c(g)+r?)(g))VPA(g) 

In lowest order perturbation theory where dt[{c{q) + rl^\q))'^q'^PA{q)]~^ = ^[(1 + 
r^'^''(g))^g^P(g)]^^ one recovers the standard perturbative result 
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provided limg2^o + ^k^\(l))'^ = 0- 

To summarize this section, we have exphcitly derived a nonperturbative flow 
equation for the scale dependence of the quark propagator ( ^.l?] ). We have also 
given a simplifled non-perturbative approximation (|5.20| ) to this evolution equa- 
tion and checked explicitly the consistency with perturbation theory in the region 
where the gauge coupling is small. In the approximation ( p. 201 ) the flow equa- 
tion ressembles a differential form of the Schwinger- Dyson equation for the full 
momentum- dependent quark propagator. In contrast to the usual Schwinger- Dyson 
equation we notice that only a small range in momenta ~ k"^ contributes to the 
momentum integral on the r.h.s. of eq. ( ^.201 ). This is therefore much easier to 
control. The actual complicated scale and momentum dependence of Z^c{q) , which 
reflects that the physical picture varies from high to low scales, arises then as a prop- 
erty of the solution of the differential equation ( p. 20 ) rather than through the form 
of the equation itself. This solution will provide the anomalous dimension rj^ and 
the function c(g) needed for the flow equation of the quark four-point function. We 
have explicitly computed the contributions 7a — e -|- 7c to the non-perturbative flow 
equation for the quark four-point function ( ^.3|) , (|5.7| ), It is straightforward to 

extract the contributions from 7^^ and 7^ from and (5.16). (The latter may 

be generalized to include Xp etc..) The result gives the non-perturbative evolution 
equation for the momentum- dependent four-quark coupling. 



6 Conclusions and discussion 

In this paper we have developed a formalism for integrating out the gluon flelds in 
order to obtain an effective action for the quarks. This is not done at once since such 
an approach would lead to complicated nonlocalities and a reliable direct computa- 
tion seems almost impossible. Instead, we account for the gluon contributions to an 
exact nonperturbative flow equation. At every scale k this needs only information 
about the "classical solution" for the gauge fleld in presence of fermions and about 
quadratic gauge fleld fluctuations around this solution with momenta k"^. The 
flow equation describes the scale dependence of an effective average action [■?/'] 
which only involves the quark flelds. In the present work we have concentrated 
on analytical work whereas the numerical exploitation of our formulae is left to a 
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separate investigator! [ITl 



Our first main result is the evolution equation (6.9) for the gluon propagator in 
the heavy quark limit. As described in sect. 3, its quantitative solution is connected 
with the heavy quark potential. The second result concerns the effective action 
for light quarks. The evolution equation for the two- and four-point functions can 
be extracted from sect. 5. The solution for the four-point function is expected 
to develop pole-like structures connected to mesons. This can be treated with the 
composite-field methods developed in ref. 0], such that one can finally make a 
transition to an effective theory for mesons describing the low momentum behaviour 
of QCD. The main difference of the present formalism as compared to ref. P] 
concerns the treatment of the gluons: While ref. accounts for the effects of 
gluons only by a phenomenologically motivated four-quark interaction, no such term 
is introduced here by hand. At short distances we simply start with the QCD action 
for quarks and gluons. The effective four-quark interaction should arise as a property 
of the solution of the flow equation. For this purpose it is crucial that the gluon 
contributions to the evolution equation are properly taken into account. 

In order to obtain the fc-dependent classical solution for the gluon as a func- 
tional of the quark fields ip we need knowledge about the effective action [■?/', A] 
for both quarks and gluons. The same holds true for the quadratic fluctuations of 
the gluon field around this solution. More concretely, the exact evolution equations 
for the fermionic two- and four-point functions involve the gluon propagator for 
t/) = 0, the ipipA, (?/'?/')^74, {tljilj)A'^, (iptp)'^ and tptpA^ vertices as well as the gluonic 
vertices and A'^ in F^ [■?/', A]. Obviously, these quantities can be computed only 
approximately and truncations are needed. We propose to use the nonperturbative 
evolution equation for F^ [■?/', A] in order to compute the inverse propagator (~ A^) 
and at least one vertex (~ A'^). (Other vertices can then be related to the A'^ ver- 
tex). Now the reader may ask why we do not work entirely in the framework of 
the evolution equation for Tk[ip, A], extracting T[iIj] only at the end of the evolution 
for k = 0. Indeed, we have recovered the perturbative /5-functions as limiting cases 
of our nonperturbative flow equations for small gauge couplings, and these pertur- 
bative /3-functions are certainly easier obtained in the framework of the evolution 
equation for Ffc[?/',74]. Also, the heavy quark potential only involves a computation 
of the gluon propagator encoded in Tk[ip, A]. The main virtue of our approach con- 
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cerns the nonperturbative aspects of the flow equations with hght quarks: We want 
to obtain a quantitatively reliable flow equation for the full momentum dependence 
of the quark four-point function. This encodes the formation of meson-bound states 
as pole-like structures in the s-channel. In turn, this requires a control of the mo- 
mentum dependence on the r.h.s. of the flow equation. Using a flow equation for 
rfe[V',A] we would need a computation of the momentum dependence of the effec- 
tive vertices tjjtjjA, {tptpyA^, (■?/'?/') A^, {tptpyA^, tjjtfjA^, and A^ in addition to the 
momentum dependence of the gluon and quark propagator and the {^ipY vertex 
in Ffc ['?/', A]. Using the effective average action for quarks [■?/'] a large part of this 
momentum dependence is encoded in the two- and four-point functions in r^f"?/^]. 
(Note that the effective {ipipY vertex in TfJ^ip] is different from the corresponding 
one in [■?/', A\ since effects of mixed quark-gluon vertices are included through the 
classical solution for A.) One may therefore hope that the momentum dependence 
of the propagator and four-quark interaction in Tk^^p] includes the dominant effects, 
whereas a less precise estimate is sufficient for the vertices appearing in the contri- 
butions from the gluon fluctuations around the classical solution. In the truncation 
used in this paper we neglect for a computation of the classical solution and the 
gluon fluctuations the vertices ipipA^, {ipipYA, {ipip^A^ and ipipA'^ and we describe 
the three vertices ipipA, A'^ and A'^ by one common scale-dependent, but momentum- 
independent, coupling constant g{k). Clearly, establishing the flow equations for the 
momentum dependence of the ipipA, A^ and vertices and using the appropriate 
solution on the r.h.s. of the flow equations should be one of the next steps in our 
approach. 

One may suspect that even the non-perturbative treatment of this paper breaks 
down for scales k of the order or below the confinement scale. Fortunately the 
formation of meson-bound states is expected at a scale considerably higher than 
the confinement scale. For the pseudoscalar mesons a first computation indicates 
k^ ^ 650 MeV In view of this one may hope that an understanding of the 
formation of mesons does not necessitate a very detailed understanding of the physics 
near the confinement scale. There are still several "hopes" and "expectations". A 
quantitative computation of the quark condensate < ipip > and the pion decay 
constant along similar lines as in ref. 0, but using the fiow equations proposed 
in the present paper should decide whether we are on a reasonable track for an 
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analytical understanding of QCD. 



Note added: Results for the form of the heavy quark potential have recently been 



obtained using flow equations [T^ and are closely related to sect. 4 of the present 
work. 



Appendix 

For the derivation of general identities it is convenient to work with a generalized 
field X 

It is composed of real gauge fields = A^^ — A^, ghosts and antighosts as well 
as complex spinors tpm and complex scalars in some representations of the gauge 
group. We use here a notation which can be employed both in coordinate space 
a{x) = a*{x) and in momentum space where a*{p) = a{—p) and = 
{2ttY6{p — p')6°'^ . For real scalars or Majorana spinors one should omit the doubling 
of arguments in x and impose f*'{p) = ip'{—p) or similar for Majorana spinors. In 
our notation the quadratic part of the action reads 

^'^ll -^^^U^P^P^^P^P^^^^P'^ = Is^SxfsXa (A.2) 

In the last part of (A.2) we have combined internal indices a and momentum labels 
p to a collective index a. The infrared cutoff is introduced as a quadratic block 
diagonal piece in the action 

with TZ^"^^* {p,p') = TZ^^\—p, —p'). We also introduce sources 

J=iK*,-C,C,-r],V,J*,j) (A.5) 
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and define 



S — So + Sgi + 5*211 + AjfeS* + S'source (A- 6) 

Here 5*0 is the classical gauge-invariant action, Sgf and Sg^ are gauge-fixing and 
ghost terms and 

For the Legendre transform = J^o'a — Wk the following identities hold 



(9F - (9F 

' Maf,Jp, T^=M^pJp, (A.9) 



a2-n 



d^Wk d^Tk d^Wk d^T 



dJadJfs dapda^ dJadJjs dapda^ ^"'^ (A. 10) 

The appearance of the matrix M = diag{l, —1, —1, —1, —1, 1, 1) reflects the anti- 
commuting properties of the Grassmann variables ^, ^, ip and ^/^ in a notation where 

o- = (a, -?/^,?/^,v2,V5*) 

a = {a*,l^,i,,^P,ip*,if) (A.ll) 

and a = A — A. Taking a derivative of (A. 6) with respect to t = In A; and noting that 
only AfcS" depends on k, one finds the flow equation for F^ = F^ — ^TZa/sC'isd-a 

r._ /^(2) _ x-i 



dtT, = -STTi^dtnk[Tl'> + nk) I (A.12) 

with STiA = TiMA, STiAB = STiBA, and Tr = / ^ Ea- The inverse propa- 
gator 

r(2) = f (2) - 7^, is given by 

We next turn to the anomalous Ward-Takahashi or Slavnov- Taylor identities for 
which we follow closely the treatment of ref. |0. For a general gauge- fixing G^{a) 
linear in (and possibly depending in addition on the background field 

Sg^ = - j rf"x?^(D^(a + A))y^c^ (A.14) 
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the sum So + Sgi + Sgh is invariant under the BRS variation x ^ X + ^brsX^j X 

X + ^BRSX^ 



^(D,ia + A)cr,-lry^cyd", -—i 

,9 2 ag 



ag 2' 



-2c^(T;)^„7/>;, -^c^(T,)^„<, -ic\T:)abipl\ic'{T:),,ip',^ (A.15) 

Here f^^^ are the structure constants of the gauge group and the hermitean 
generators in the appropriate representations. It is useful to introduce external 
sources = 

P = (^^,-7^7^-5W,5W^^M,5M) (A.16) 



for the BRS variations of x or X ^ind to define VFfc[J, /?] similar to (|A.6| ) by adding 
in ([A .71) an additional source term 

SSlrce = -/5a(5BRSX)a = -/^^(^BRSX)^ (A. 17) 

The BRS-invariance of the measure implies 

= < f^BRS'S' > = < 'R-al3Xl3{.^-SKSX)a >[7=0 " < -'a('^BRSX)a >|7=0 

= < 'Tlapi^BKSX) pX-yM^ta >|7=0 " < Ja{5BKSX)a >|7=0 (A. 18) 



or 



Ol3a 1^=0 \OJl3 OJp ) dl3a |^=o 



Using ( |A.8| ), ( |A.[j| ) and the identities 



(A.20) 
(A.21) 
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one arrives at the identity 

^MM^^^^U)-!^^^^^^ (A.22, 

In addition, we note the simple identity 

^ = —G\a) (A.23) 

which allows to eliminate the source 7. Similarly, linearity in c yields the field 
equation for the antighost 

where 

G7=d^ = iDM)yy (A.25) 



depends only on the background field A. For the last identity in ( [A .251) we have 



used the particular gauge-fixing = [D^(y4)]^^a^, which will be assumed in the 
following. For this gauge we now insert (A.23), (A. 24) such that the Ward identity 
reads in explicit components in momentum space P] 

d'^p ( dV dV dV dV 



d^m (P) dES'^ ip) d^Prn (p) d& (p) 

+ . = 41s + 4ts (A.26) 



with 



d^aip) d6^J^\p) d^*aip) d6i^\p) 
d'^p dy d'^q i^(A)uur /N/-nf2) 



1 f d'^p dy 



-^y(2^(4^i^^^^S(^^pO(^(^^ + 7^)-,^.,^ 
-^/-»/ (05^^(^-p)^S(P'PO(r(^) + 7^)^4,)^^(^)| (A.27) 
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and 



Here we have subtracted the "bare" gauge-fixing term 

r' = rfc|^^„ - ^ 1 d''x{D^{A){A^ - A,,nD,{A)iA, - A,)y (A.29) 

and we use the index convention (cf. ( |A.13D )p| 

(2) _ d^T (2) _ d^T 



p(_2) ^ -^(2) 



i^ipmq) ^(qWp) d^{q)dtlj{p) d^{p)dtlj{q) 
(2) ^ d^T (2) ^ 

We recover the usual identities in the hmit k ^ since TZ = imphes a vanishing 
BRS-anomaly ^brs = 0. 

The sources jS appear only linearly in 5* and it is straightforward to derive the 
identity 

dT dW dS 1 1 

= =< >= --DjAyy^y - -f^^y < a^cy (a.si) 

where the connected two-point function < ac >c is related to the appropriate matrix 
element of the propagator (F^^) + n)-^ by { \KA^ ) 

< a;cy >,= (r(2) + n)-l^, (A.32) 

This allows one to eliminate the explicit dependence on the source (5 in favour of 
expressions containing two- and three-point functions and to restrict the discussion 



^Note the minus sign whenever the second index is ^ or tp. 
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to /3 = afterwards. Identities similar to ( |A.31| ) for the other sources of BRS- 
variations are easily derived 

dr 

OOa 

i{T:u{e^i+ < e^'i >) (A.33) 



Finally, one has the field equation for the ghost field 

+D^{AYyD,{A)y'^C - < alD^iA)y^c'" (A.34) 

For an abelian gauge theory the vanishing structure constants f^yw lead to im- 
portant simplifications: First one can replace (|A.31|) 

^--W, f-O (A.35) 

di3t, g <97 

and evaluate everything for 7 = = 0. The field equations for the ghosts reduce 
to 

dr 
dr 

where we have limited the matter content to a complex scalar with {Tz)ab = —^ab- 
The ghost-dependent part of F is therefore uniquely determined 

Fgh = / d^x [d^ld^i - iC5ip + 1^5^*} (A.37) 

It is independent of k and equals the classical expression. Eq. ( [A.37| ) also implies 
that for 5 = 5 = the connected two-point functions involving ^ and if' , ip*' vanish 

< C ^' >C|,^o = < c >e|,^„= (A.38) 

and therefore 

00 \S=S=0 00 \5=5=0 



-d.d^i (A.36) 
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We can now evaluate the Ward identiy ( |A.26| ) for = 0, 6 = 6 = where A^^^ = 0: 



dV 



+ 



-(p*{p + q) 



d(p*{p) 

^dpl 



cfip 
(27r)« 

_ /iM 



dip{p) 



^{p - q) 



(A.40) 



BRS 



{7^(-)*b,p')(^(^) + 7^);(;,) 



{2nY {2nY (27r)'^ ^ ' ^ v^' ^ / ' ' -/^(p')^(P-g) 

-7^(-)(p,p')(^(^) + 7^)^l^,)^(^^^)}e(g) (a.4i) 

Since for 6 = 6 = the ghost sector decouples completely from the {ip, (f*, A) sector 
we can evaluate (F^^^ + TZ)~^ for ^ = ^ = 0. Furthermore the definitions ( [A.l(]| ), 
( lATll) imply 

(r^^^ + =< ^*'ip')v\p - q) >c= (r(^) + 7^)^J,_,),(,,) (A.42) 

and n^^>{p,p') = n^'P\p',p). Since ( [A .401) must hold for arbitrary ^(g) we finally 
obtain the modified Ward identity in a form not involving the ghost anymore 



f'dA,{q) ' J i2n 



+ 



d'^p 



dV 

dip{p) 



^{p - q) 



dip*{p) 



V*{p + q) 



(A.43) 



d'^p dy 



[{n^^\p,p' + g) - n^^\p - g, J9'))(F(2) + 7^)-(;,) 



(27r)^ (27r)'^ L^'^ ^■^'■^ ' ■^'/-//v- ' '-J^(p'Mp)\ 

The existence of such a form is not surprising since for our choice of gauge fixing the 
ghost sector is just a free field theory and can be omitted altogether. We observe 
that the l.h.s. of eq. (|A.43| ) is simply a local gauge variation of F' with A kept fixed. 
The gauge-invariant and A- independent part of F' does therefore not contribute. 
From Tk{f, A, A) = Tk{ip, A, A,^ = ^ = 0) we may subtract a gauge-invariant 
kernel Tjs[ip, A] = Tk[if, A, A = A] and define 

fgf = Ffc A,A]-Tk[f,A]- Fgf [A, A] (A.44) 

We can therefore replace F' by Fgf on the l.h.s. of (|A.43|) . The Ward identity only 
constrains the "generalized gauge-fixing term" Fgf ^ which contains the /c-dependent 
counterterms for k > and vanishes for k = 0. The invariance of the average action 
with respect to simultaneous gauge transformations of A and A implies 



1 



or 



+ 



or 



dA,{q) dA,{q)^ 



+ 



d'^p 
(2^ 



^F 

dip{p) 



^{p-q) - 



or 



dip*{p) 



^*{p + q) 



= 
(A.45) 
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and similarly for V or f gf leading to an alternative form of the identity ( |A.43| ) 

(A.46) 

It is interesting to compare this equation to an identity for the background field 
dependence of f gf derived earlier 

dMqj ~ J {27rY{27rY dA,{q) + ^hiP'Mp) ^^"^^^ 

One can show that the Ward identity |A.46| ) can be derived from the more 
general "background field identity" (|A.47| ) by using the invariance of A]^^ S under 
simultaneous gauge transformations of (f, A and A which yields 

I dV^pP^ = ^^^Hp - q,P') - 7^(^HP,P' + q) (A.48) 
9 dA^il) 



The background field identity ( |A.47|) has a simple solution in the approximation 



where the A dependence of the propagator (F*^^) +TV) ^ on the r.h.s. is neglected 

fgf =< a'^^S > -Al^^SM = / ^d-0j-X"^(P^P') < ^*iPMp') >c (A.49) 
A background field identity for non-abelian gauge theories has also be derived |H, 



TT]] . The precise relation to the Slavnov- Taylor identity has not yet been established. 
It is clear, however, that the background field identity contains information beyond 
the Slavnov- Taylor identity. 

Let us finally turn to solutions of the Slavnov- Taylor identity for non-abelian 
gauge theories. We concentrate first on a vanishing infrared cutoff TZk = 0. We are 
interested in solutions for vanishing sources f3 for the BRS- variations. Eq. (A. 26) 
contains then a sum of expressions for which the first factor can be evaluated at 
f3 = 0, whereas the second factor involves the coefficients linear in f3 evaluated at 
j3 = 0. It is straightforward to show that the ansatz 

r = r,„4AV',<^] + rgf + rgh + r, + Ar[A] 

Tgf = ^ld'x{D^{A)iA,-A,)nD'^iA)iA^-A^r 
Tgh = - J d'xe{D^{A)D^{A)0' 
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+^Si^^C{T.)am - ^Si^k%T:Uvl} (A.50) 

obeys the identity (A. 26) with vanishing anomaly on the r.h.s. Here rinv[A,ip,(f] is 
an arbitrary gauge-invariant functional which does not depend on the ghost fields, 
and Ar[y4] is an arbitrary gauge- invariant functional of the background field A. In 
fact, gauge invariance of imphes the identity 

(T.).,^,)(p) - ^^^AeiT:U^l)ip)} = (A.51) 

This allows to replace F' by Fgh for the first factors in (A. 26). What remains is the 
relation 

[D^{A)0'{p) - ^r'-^^mnip) = o (A.52) 



/ 



which is easily verified using the Jacobi identity for the structure constants f^*^ p'^v — 
jzty jtsw _ jzstjtyw ^ obscrvc that the gauge- invariant part Fj„t, remains com- 
pletely unconstrained by the Slavnov- Taylor identity. We can also verify that the 
ansatz ( [A. 50] ) obeys the field equation for the antighost (A. 24). In contrast, the 
source identities (A. 31), (A. 33) hold only if the pieces involving the connected two- 
point functions < ac >c, < cc >c, < ap' >c etc. all vanish.^ This is generically not 
the case for the ansatz ( [A.50| ), since the cubic vertex ~ g^S,A in Fgh induces a non- 
trivial off-diagonal matrix element in the inverse propagator F^ , mixing the ghost 
sector to the other fields. This, in turn, is responsible for corresponding off-diagonal 
elements in (F*^^^ -|- 1Z)~^. We observe that these off-diagonal elements vanish in 
the limit of small gauge coupling (? — such that (A. 31) and (A. 33) are obeyed in 
this approximation. The situation for the ghost field equation (A. 34) is completely 
analogous. We conclude that for small g ^ the whole picture becomes formally 
very similar to the abelian case discussed above, with leading non-abelian structure 

^These two-point functions should be evaluated from the ansatz ( A.50|) for vanishing sources [3. 
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given by the ansatz ( [A. 501 ). We should emphasize that this approximation can be 
expected to be vahd only for the momentum range where the running renormalized 
gauge coupling remains small enough. 

The ansatz ( |A.50| ) is, however, not the most general solution. We may try 



to obtain a more general solution by replacing in (|A.5CI|) the ghost field C,^ by a 



functional ^^[A,^,^,ilj,ip] and changing the term linear in 7^. The functional 
should have ghost number one and obey the same symmetry transformation laws as 
under gauge transformations acting on A and A simultaneously. Eq. ( [A.51| ) holds 



also with ^ replaced by ^ such that Tinv again drops out. The remaining equation 
relates |^ to the functional form of ^ 



2 goA'li^p) 



d^mip) dip^ip) 

difiaip) 



-^^Mii'iTy)^m){p) + '^^^(^'(T;Uvl)ip)} = (A.53) 



Given an arbitrary ^ one can always solve the equation for |^ provided the rela- 
tion between ^ and ^ remains invertible. At this level we have therefore a general 
class of solutions for the identity (A. 26) since ^ remains essentially unconstrained. 
The constraints arise from the source identities (A. 31), (A. 33) and the ghost field 
equation (A. 34). (Note that the field equation for the antighost (A. 24) is obeyed 
for arbitrary In fact, the difference between ^ and ^ is related to the connected 
two-point functions involving the ghost field 

(T,)„5(r - Dv^fc- < c^;, >c) = (A.54) 

Similar relations must hold for v^*,^' ^-^id and the equation for |^ must be com- 
patible with the solution of ( |A.53| ). We conclude that the generalized ansatz with ^ 
replaced by ^ can only be used in the approximation where this system of equations 
for ^ is self-consistent. Beyond this approximation the gauge invariance of the sector 
with ghost number zero {^mv^ can probably not be maintained. 
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In the presence of a nonvanishing infrared cutoff TZk an additional piece 
containing counterterms should be added. It vanishes in the limit k ^ 0. In the 
perturbative regime contains a gluon mass term ~ fc^ (see sect. 4). There is 
no reason why Td should be gauge-invariant under transformations which leave the 
background field A fixed. We may define Td by the requirement that F — Td obeys 
the anomaly-free Slavnov- Taylor identity. The anomaly ^brs + ^brs '^'^ the r.h.s. 
of eq. (A. 26) determines then the form of Tct- In a lowest order approximation we 
may treat Td and the anomaly as a small quantity and linearize eq. (A. 26) in Td- 

For the computation of the gluon propagator in the present paper we make only 
a crude approximation for the three- and four-gluon vertices appearing on the r.h.s. 



of the flow equation. They correspond to the ansatz (|A.5CI|) with Fj„^ containing 
only a piece ~ F^^F^" . An improved treatment of these vertices could generalize 
Tinv for a nontrivial momentum dependence of the propagator 

Tinv = \j dxF^,K{-D\A))F^'^ (A.55) 

with K{x) = {Ga{x) - Ga{0)))/x. 
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Figure Captions 

Fig. 1: Wave function renormalization Zp as a function of the average scale k. 

Fig. 2: Momentum-dependent anomalous dimension x(g) in comparison with per- 
turbation theory for various values of A; > 1 GeV. 

Fig. 3: The same as fig. 2, for /c < 1 GeV. 
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